A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 

06(02.01.2; 08.02.1; 08.09.2; 08.14.1; 13.25.5) 



ASTRONOMY 

AND 
ASTROPHYSICS 



On the X-ray fast-time variability of Sco X-2 (GX 349+2) 



o 
o 

(N 

X) 
CD 



> 
O 



Oh' 
I 

O 



Paul M. 0'Neill\ Erik Kuulkers^'S, Ravi K. Sood\ and Tadayasu Dotani'* 

^ School of Physics, Australian Defence Force Academy, University of New South Wales, Canberra ACT 2600, Australia 

(p.m.oneill@adfa.edu.au r.sood@adfa.edu.au) 
^ Space Research Organization Netherlands, Sorbonnelaan 2, 3584 CA Utrecht, The Netherlands (erikk@sron.nl) 
^ Astronomical Institute, Utrecht University, P.O. Box 80000, 3507 TA Utrecht, The Netherlands 

Institute of Space and Astronautical Science, 1-1 Yoshinodai 3-chome, Sagamihara-shi, Kanagawa 229-8510, 

Japan (dotani@astro.isas.ac.jp) 



Received / Accepted 



Abstract. We have analysed archived Ginga data on the 
Z source Sco X-2 (GX 349-^2). We present the first de- 
tailed investigation of its X-ray fast-time variability, as a 
function of position in the Z track. During the two-day ob- 
servation over the period 5-7 March 1989, the source was 
in the so-called flaring branch, and the lower part of the so- 
called normal branch. We found broad peaked noise with 
a centroid frequency and width of ~4-7 Hz and ~6-12 Hz 
respectively. The peaked noise was strongest in the lower 
flaring branch, with a maximum fractional rms amplitude 
of ~3 %. We conclude that it is not a manifestation of atoll 
source high frequency noise, as had been suggested, and 
compare it with the power spectral features seen in other Z 
sources. We find that the peaked noise is markedly differ- 
ent to the quasi-periodic oscillations found in the normal 
and flaring branches of Sco X-1; however it bears some 
resemblance to that seen in the flaring branch of Cyg X-2 
at low overall intensities. 

Key words: accretion, accretion disks - stars: binaries: 
close - stars: individual: Sco X-2; GX 349-1-2 - stars: neu- 
tron - X-rays:stars 



1. Introduction 

Sco X-2 (GX 349+2) is a neutron star low-mass X-ray 
binary and Z source. Z sources trace out a characteris- 
tic 'Z' shaped track in an X-ray colour-colour diagram 
(CD) or hardness-intensity diagram (HID) consisting of a 
'horizontal' branch (HB), 'normal' branch (NB) and 'flar- 
ing' branch (FB) (Hasinger & van der Klis 1989; hereafter 
HK; van der Klis 19951 ). The orientation of each branch in 
a GD or HID depends upon the choice of energy bands. 
Mass-accretion rate (M) is thought to increase from the 
HB, down the NB and up the FB. The other traditional Z 
sources are: Sco X-1, GX 17+2, GX 340+0, Cyg X-2, and 
GX 5-1. 



A class of objects that are related to Z sources are the 
'atoir sources. Atoll sources trace out a curved track in 
a CD, with inferred mass-accretion rate increasing along 
the track. Atoll sources are thought to have lower mag- 
netic field strengths, and lower mass-accretion rates than 
Z sources (HK). 

The X-ray fast timing properties of Z and atoll sources 
are different, yet related. The features observed in the 
power spectra of both classes are normally correlated with 
the position of the source in a CD or HID. Z sources 
exhibit very low frequency noise (VLFN), low frequency 
noise (LFN), high frequency noise (HFN), and quasi- 
periodic oscillations (QPOs). QPOs with frequencies in 
the range 15-60 Hz (called HBOs) and 200-1200 Hz (kHz 
QPOs; these typically occur in pairs) are present in the 
HB and upper part of the NB. Atoll sources exhibit VLFN, 
HFN, 1-70 Hz QPOs and kHz QPOs. AtoU source HFN 
can sometimes be peaked, and resemble broad QPOs (HK 
and references therein). A recent review of the fast-time 
variability of Z and atoll sources, with a discussion of var- 
ious models, has been given by van der Klis (200C). 
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In Z sources, narrow QPOs (called NBOs) with cen- 
troid frequencies in the range 5-8 Hz have been observed 
in the NB. NBOs typically become detectable halfway 
down the NB, and their properties differ between ob- 
jects. In Sco X-1 the frequency of the NBO is ^5.5-8 Hz 
in the NB, and jumps to ^^12 Hz as the source enters 
the FB (the QPO is th en referred to as an FBO) (Di- 
eters & van der Klis 200C). As the source moves up the FB, 
the frequency and full width at half maximum (FWHM) 
increase; eventually, at ~10 % of the way up the FB, the 
FBO becomes so broad as to be indistinguishable from 
the underlying noise (but the relative width never exceeds 
~0.5). In Sco X-1 the fractional rms increases with in- 
creasing M and is greater at higher energies; and there 
is no time lag b etween low and high energy bands (Di- 
eters et al. ^m% . The N/FBO seen in GX 17+ 2 beh aves 
in a similar fashion to Sco X-1 (Penninx et at. 1990| ). In 
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the other Z sources NBOs have also been observed but are 2. Observations 



generally less prominent (e.g. Kuulkers et al. 1994). 

Sco X-2 is seemingly very similar to Sco X-1. The Z 
tracks have nearly identical shapes (except for the absence 
of an HB in Sco X-2) and the lightcurves of both objects 
exhibit strong flaring behaviour (Schultz et al. 1989; Ku- 
ulkers & van der Klis 1995a). The orbital periods are also 
similar: 18.9 h for Sco X-1 (Co wley fc Crampton pJ75| ); 
and ~22 h for Sco X-2 (Wachter |1997D. One might there- 



fore expect the X-ray fast-time variability in Sco X-2 to 
resemble that seen in Sco X-1; but this is not the case. 

EXOSA T ob servations of Sco X-2 reported by Pon- 
man et al. (1988) showed a complete FB. No narrow NBO 
was found; however broad peaked noise was detected with 
a centroid frequency and FWHM of about 5 Hz and 10 Hz 
respectively. The peaked noise was strongest in the lower 
FB, with a fractional rms of 3.2±0.2 %. There was no 
significant change in the frequency or width of the peaked 
noise, and there was no time lag between the low and high 
energy bands. We note that the authors divided their data 
into four segments on the basis of intensity, and not ac- 
cording to position in the Z track. By using that method 
they could not differentiate between the NB and the lower 
FB. 

A small set of Rossi X-ray Timing Explorer observa- 



tions were reported by Kuulkers & van der Klis (1998; 
hereafter KK), when Sco X-2 was near the NB/FB ver- 
tex. In the lower NB they found peaked noise with a cen- 
troid frequency, FWHM and fractional rms of 9.4±0.5 Hz, 
16±2 Hz and 2.7±0.1 % respectively. In the lower re- 
gion of the FB they measured these values as 5.8±0.2 Hz, 
11.0±0.7 Hz and 4.2±0.1 %. There was a small but signifi- 
cant decrease in frequency when the source moved through 
the NB/FB vertex. They commented that the peaked 
noise resembled the peaked HFN sometimes seen in atoll 
sources. They found that the peaked noise fractional rms 
was stronger at higher energies. 

Even though there have been clear indications that 
Sco X-2 exhibits marked differences to Sco X-1, it has 
remained largely ignored. This is perhaps due to the fact 
that it looks so similar to Sco X-1 in a CD and HID. 

We present results obtained from an analysis of 
archived data from the Ginga satellite; these data have not 
previously been presented. This work is the first detailed 
analysis of the power spectral properties of Sco X-2 as a 
function of position in the Z track. This analysis method 
is extremely useful because it allows us to unambiguously 
distinguish between the NB and lower FB, and is help- 
ful when making comparisons between different objects. 
We further highlight the differences between Sco X-2 and 
Sco X-1, and conclude that the observed peaked noise is 
not the same as atoll source peaked HFN, as had been 
suggested. We compare our results with the power spec- 
tral features seen in other Z sources. 



The Ginga Large Area Counter ( LAC) 
(Makino & ASTRO-C Team |1987| ; Turner et al. [1989^ 
observed Sco X-2 between 5 and 7 March 1989. We re- 
trieved the data from the Leicester Database and Archive 
ServiceQ. In Table |l| we present a log of observations. 

The LAC consisted of eight detectors, and was sen- 
sitive to X-rays in the range 1-37 keV. Several modes 
were available, with differing time and energy resolutions 
(higher time resolution sacrifices energy resolution). The 
observations discussed here were done in the so-called 
MPC3 and PC modes. For the MPC3 mode the coarse 
gain setting was high, providing an energy range 1-19 keV 
over twelve channels. The time resolution was 7.8 ms and 
62.5 ms for high (MPC3-H) and medium (MPC3-M) bit 
rates respectively. In the PC mode, the LAC was divided 
into two groups of four detectors each. The coarse gain set- 
ting was low, giving an energy range 1.5-24.3 keV for one 
detector group, and 1.2-17.9 keV for the other. Each de- 
tector group had two energy channels, so the PC mode ob- 
servations had four overlapping energy bands. In the high 
bit rate mode (PC-H), the time resolutions were 0.98 ms 
and 1.95 ms for the lower and higher energy bands re- 
spectively. Our 0.98 ms data were rebinned with a time 
resolution of 1.95 ms, to obtain lightcurves for the full 
energy range. 

During the observations in section C (see Table P the 
high voltages in detectors 5 and 6 were accidently changed. 
The energy range for the MPC3 observations stated above 
is therefore only approximate for section C. 



3. Analysis 

The counting rates were corrected for detector deadtime 
and collimator response. During the observations in sec- 
tions A-C, the collimator response was ~0.96. In section 
D the response was ~0.73. The data were not corrected 
for collimator reflection. The effect of collimator reflection 
was to overestimate co untin g rates below ~6 keV by a few 
percent (Turner et al . 1989; see also Hasinger et al 
|1997l ). 



1990 



Wijnands et al. 

Hardness-intensity diagrams were constructed using 
64 s averages. For the MPC3 data we defined the hardness 
as the counting rate ratio between 8.2-16.7 keV and 5.8- 
8.2 keV. The intensity was defined as the counting rate in 
the range 1.1-16.7 keV. For the PC data we defined the 
hardness as the counting rate ratio between 5.7-17.9 keV 
and 1.2-5.7 keV. The intensity was defined as the counting 
rate, from one detector group only (i.e. four detectors), in 
the range 1.2-17.9 keV. 

Timing analysis was performed on the raw (i.e. not cor- 
rected for deadtime and collimator response) data. Power 
spectra were calculated from 16 s intervals in several en- 
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Table 1. Observation Log 



Section 


Start Time (UT) 


End Time (UT) 


Usable Data (s) 


Modes 


Time Res (ms) 


A 


5 Mar 1989 20:58 


6 Mar 1989 00:17 


4512 


PC-H 


0.98/1.95 


B 


6 Mar 1989 00:27 


6 Mar 1989 19:13 


2880 


MPC3-H 


7.8 








12928 


MPC3-M 


62.5 


C 


7 Mar 1989 02:13 


7 Mar 1989 14:30 


12928 


MPC3-M 


62.5 


D 


7 Mar 1989 18:12 


7 Mar 1989 19:14 


1184 


PC-H 


0.98/1.95 



ergy bands, and then averaged according to either cor- 
rected counting rate or rank number. 



(see also Lewin et al. 1992; Hertz et al 



et al. 


1990 


1992; 


Di- 



eters & van der Klis pOO0| ) and is a one-dimensional mea- 



sure of the position of the source in a CD or HID. Rank 
number increases smoothly with increasing inferred mass- 
accretion rate. To measure the rank numbers, a spline was 
fitted through the Z track, and each point was projected 
onto it. Rank number was then determined by measuring 
the distance of each point along the spline. We defined 
the NB/FB vertex as rank number 1, and the top of the 
observed FB as rank number 2. 

We calculated the expected deadtime-affected white- 
noise levels and then subtracted the white-noise from 
each average power spectrum. The procedure to calcu- 
late white-noise levels in sum channel power spectra is 



well known, and is described by van der Klis (1989) and 
Mitsuda & Dotani (1989). The full procedure to calculate 



white-noise levels in energy resolved power spectra is not 
explicity stated in those references; we therefore describe 
our procedure in the Appendix. 

After the deadtime-affected white-noise levels had 
been subtracted, the spectra were normalised to 
fractional-rms-squared per Hz and functional fits were 
performed. Very low frequency noise (VLFN) was mod- 
elled using a power law A v°' and integrated over the range 
0.01-1 Hz. Peaked noise was modelled using a Lorentzian. 
All power spectra were initially fitted with only a power 
law. The significance of adding a Lorentzian component 
to the fit was then estimated via an F-test. If the in- 
clusion of a peaked noise component was not significant 
(< 2>a) we measured 90 % confidence upper limits by using 
= 2.71. Errors on all parameters were measured us- 
ing Ax^ — 1. When a fit was not sensitive to the centroid 
or FWHM of the Lorentzian (or was not significantly im- 
proved by its inclusion) we fixed those parameters at val- 
ues measured from other power spectra (the actual values 
are shown in Table ^). We fitted each PC mode power 
spectrum both with and without an extra free constant 
component. The reasons for doing this are discussed in 
Sect. 



4.2 



All fractional rms amplitudes were corrected for chan- 
nel cross-talk (Lewin et al. 1992). We assumed varia- 



tions at all energies were in phase. Finally, the fractional 
rms values w ere m ultiplied by a binning correction factor 
(van der Klis |l98g| ). 



5000 



4 .- * 



10' 

n tensity (counts/s/8- detectors) 



1 .5x1 0' 



Fig. 1. HID of section B, with 64s averages. The intensity 
is the counting rate in the range 1.1-16.7 keV. The hard- 
ness is the counting rate ratio between 8.2-16.7 keV and 
5.8-8.2 keV. The error bars are slightly larger than the 
dots. 



4. Results 

Hardness-intensity diagrams 

Fig. |l| shows a HID for all section B data. The small, 
approximately vertical NB is clearly visible. Sco X-2 is 
known to flare to intensities roughly twice the level of the 
persistent emission. The highest intensity in Fig. |l] is a 
factor of 2.2 greater than the lowest intensity: we conclude 
therefore, that we have observed a fully developed flaring 
branch. The results from section C are compatible with 
those from section B; however due to the uncertain energy 
ranges of those data, we will not discuss them any further. 

A HID of sections A and D is shown in Fig. ||. The 
lowest and highest intensities in Fig. || differ by a factor 
of 2; therefore (in comparison with Fig. |^), nearly all of 
the FB was traced out. No clear NB/FB vertex is ap- 
parent (but we were able to estimate where it might lie; 
see Sect. When we constructed a HID from MPC3 
mode data using energy bands similar to those available 
in PC mode, the presence of the NB was much less pro- 
nounced; the NB became more like a 'blob', rather than 
a clear branch. Therefore we caution that in Fig. |^ we 
cannot confidently distinguish between data points from 
the lower NB (if there are any) and those actually in the 
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Fig. 2. HID of sections A and D, with 64 s averages. The 
intensity is the counting rate in the range 1.2-17.9 keV. 
The hardness is the counting rate ratio between 5.7- 
17.9 keV and 1.2-5.7 keV. The error bars are slightly larger 
than the dots. 



lower FB. However, as will be seen in Sect. 4.3, our main 
interest (with the PC data) is with the top half of the FB, 
so this ambiguity is not a major problem. 

4-. 2. Timing analysis of section B MFCS mode data 

During the observations in MPC3-H mode, the source was 
found only in the lower FB. We divided the data into two 
segments with mean ranks 1.09±0.02 and 1.26±0.05. An 
example of a power spectrum (1.1-16.7 keV) is shown in 
Fig. ^, with mean rank number 1.09. The results of fitting 
the MPC3-H data (1.1-16.7 keV) are shown in Table |. 
The FWHM of the best-fit Lorentzian is greater than half 
its centroid frequency, so the feature cannot be formally 
described as a QPO. The MPC3-H data show that there is 
no significant change in the centroid frequency or FWHM 
of the peaked noise, as the source moves up the FB (the 
decrease in the centroid frequency between rank numbers 
1.09 and 1.26 is significant only at the ~2 a level). The 
fractional rms amplitude of the peaked noise is also con- 
sistent with remaining constant. 

We investigated the energy dependence of the peaked 
noise fractional rms, by calculating energy resolved power 
spectra from the MPC3-H data. We used four energy 
bands: 1.1-3.4-5.8-8.2-16.7 keV. All power spectra from 
a particular band were averaged together, giving a mean 
power spectrum from the lower FB; the mean rank num- 
ber was 1.18±0.05. We determined the centroid frequency 
and FWHM of the peaked noise feature from an average 
power spectrum from the 1.1-16.7 keV band: they were 
6.7±0.4 Hz and 8.3j;}'^ Hz respectively. The centroid and 
width were fixed at those values when fitting the power 
spectra from the other four energy bands. The energy de- 
pendence of peaked noise fractional rms is shown in Ta- 




100 



Frequency (Hz) 



Fig. 3. MPC3-H power spectrum from the energy range 
1.1-16.7 keV. The mean rank is 1.09 (lower FB). The sohd 
line shows the best fit. 

Table 3. Energy dependence of peaked noise fractional 
rms amplitude 



Energy Range (keV) 


rms (%) 


1.1-3.4 


<2.8 


3.4-5.8 


<2.6 


5.8-8.2 
8.2-16.7 


4.8lg-j 



ble |. The peaked noise is clearly stronger at higher ener- 
gies and is consistent with being constant above ~6 keV. 

The MPC3-M data from the track shown in Fig. |^ were 
divided into eight segments. Power spectra were averaged 
on the basis of rank number. An example (1.1-16.7 keV) 
is shown in Fig. ^ with mean rank number 1.09. Table |^ 
shows the results of fitting the MPC3-M power spectra 
(1.1-16.7 keV) from section B. It was only in the lower 
FB that we could meaningfully constrain the centroid and 
width of the peaked noise, which were ~4 Hz and ^6 Hz 
respectively. There is no significant change in those values 
between mean ranks 1.09 and 1.27 (Table ||). The cen- 
troid frequency is lower in the MPC3-M data than in the 
MPC3-H data. This discrepancy is probably due to a lack 
of information above 8 Hz in the MPC3-M data. 

The dependence of peaked noise and VLFN fractional 
rms on rank number is shown in Figs. || and ^ respectively. 
The peaked noise fractional rms is greatest in the lower FB 
and becomes undetectable (upper limit 1.3 %) by about 
halfway up the FB. Of note is the fact that it is actually 
detectable up to a mean rank of 1.41. (This is important 
because N/FBO typically become undetectable by ^10 % 
of the way up the FB.) The drop in power as the source 
moves further up the FB or into the NB may be due to 
the feature moving above the Nyquist frequency of these 
data; however we can exclude this effect with respect to 
movement up the FB (see Sect. O). The VLFN fractional 
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Table 2. Results from fits to the section B MPC3 power spectra from tlie energy range 1.1-16.7 keV. 



Mean Rank VLFN Peaked Noise Reduced X^/dof 
Index rms (%) Frequency (Hz) FWHM (Hz) rms (%) 



MPC3-H data 


















1.09±0.02 


1 9+0.2 
-r.Z_() 3 


1 9+0.2 


7.3±0.4 


7 1+1-2 

'•1-1.0 


2.9±0.2 


1.25 


/ 


43 


1.26±0.05 


-2 0+"'^ 


^■"-0.4 


c =+0.8 
^■•^-1.0 


10.3^^2;! 


3.0±0.3 


0.79 


/ 


43 


MPC3-M data 


















0.95±0.02 


-1.5±0.3 


1.2±0.2 


4.3=^ 


6.1^ 


1.8±0.2 


0.85 


/ 


25 


1.00±0.02 


-1.2±0.2 


1.3±0.1 


4.3^ 


6.1'' 


2.3±0.1 


1.28 


/ 


25 


1.09±0.03 


-1.5±0.3 


1.4±0.2 


4.3±0.4 


6.ii?i 


2.8±0.3 


0.94 


/ 


23 


1.27±0.06 


-1.9±0.2 


2.8±0.5 


4.4±0.6 


6.21?:^ 


2.3±0.3 


1.40 


/ 


23 


1.41±0.03 


-1.6±0.1 


2.9±0.2 


4.4=' 


6.2^ 


1.4±0.2 


0.84 


/ 


25 


1.50±0.03 


-1.5±0.1 


2.8±0.2 


4.4a 


6.2'' 


<1.3 


1.06 


/ 


26 


1.60±0.03 


-1.6±0.1 


3.3±0.2 


4.4a 


6.2'' 


<0.9 


1.10 


/ 


26 


1.75±0.09 


-1.8±0.1 


4.8±0.3 


4.4^ 


6.2^ 


<1.2 


1.10 


/ 


26 



"Fixed parameter 




0.1 



Frequency (Hz 



Fig. 4. MPC3-M power spectrum from the energy range 
1.1-16.7 keV. The mean rank is 1.09 (lower FB). The sohd 
line shows the best fit. 



rms is seen to increase as the source moves up the FB, as 
expected for Z sources (HK). The slope of the power law 
appears to be uncorrelated with rank number. 

4-.3. Timing analysis of PC mode data 

The NB/FB vertex was not unambiguously observed, so 
we did not calculate rank numbers for each individual 
interval; averaging was done on the basis of intensity. 
The track was divided into four segments; the intensity 
ranges were: 2800-3200, 3200-3900, 4500-5200 and 5200- 
5900 count/s (four detectors only). The mean intensities 
of those four segments were: 3041±6, 3458±14, 4962±25 
and 5492±31 counts/s. We could, however, estimate rank 
numbers by comparing Fig. ^ with Fig. |l| (assuming the 
intensities are derived from the same energy band). As 
stated in Sect. 4.1, the track shown in Fig. is not quite 



0.8 



.2 1 .4 

Rank Number 



Fig. 5. Peaked noise fractional rms values as a function 
of rank number for the MPC3-M data (solid dots, 1.1- 
16.7 keV), and PC data (open squares, 1.2-17.9 keV). The 
fractional rms values of the PC data are from fits which 
included a constant component. 



a complete FB. The lowest intensity observed in Fig. |^ 
is ~2900 counts/s; this can be taken as the upper limit 
of the intensity of the NB/FB vertex. Alternatively, if we 
assume the highest intensity present in Fig. || is actually 
the top of the FB, then the unseen NB/FB vertex can be 
placed at ^^2780 counts/s; this can be considered to be the 
lower limit. By considering these two extremes we calcu- 
lated the approximate rank numbers of the four segments 
as: 1.06, 1.20, 1.68, and 1.80 (±0.03). We stress that these 
rank numbers are only a rough guide, and that the lower 
FB in Fig. || may infact contain data from the NB. 

An example power spectrum is shown in Fig. ^ (1-2- 
17.9 keV), with mean intensity 3041 counts/s (mean rank 
number ~1.06). The results from fitting the PC power 
spectra (1.2-17.9 keV) are shown in Table ^, where we 
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Fig. 6. VLFN fractional rms values as a function of rank 
number for the MPC3-M data (soHd dots, 1.1-16.7 keV), 
and PC data (open squares, 1.2-17.9 keV). The fractional 
rms values of the PC data are from fits which included a 
constant component. 



1, il 



Frequency (Hz) 

Fig. 7. PC power spectrum from the energy range 1.2- 
17.9 keV. The mean intensity is 3041 counts/s (mean rank 
number ~1.06; lower FB). The solid line shows the best 
fit, and includes a constant component. 



present the fits from models both with and without a con- 
stant. The results from fits with a constant are shown 
also in Figs. || and The most important thing to note 
is the nondetection of peaked noise in the upper FB (up- 
per limit 1.6 %). The PC data Nyquist frequency is 256 Hz; 
therefore the nondetection of peaked noise in these data is 
not due only to the feature moving to higher frequencies. 
From this, we conclude that the nondetection of peaked 
noise in the MPC3-M data (Table |^) is also a real weak- 
ening of the peak. 

The reduced x^s of the fits without a constant in Ta- 
ble ^ are formally acceptable; however visual inspection 
of the power spectra revealed excess power above about 
100 Hz (e.g. see Fig. @), and the best-fit VLFN indices 
found in the two lowest intensity segments are unusu- 
ally flat. Given these indications, we attempted to account 
for a possible high frequency feature. We added a cut-off 
power law to our model and attempted to fit the spectrum 
from the lowest intensity segment. We could not constrain 
the index or cut-off frequency of this component. When 
we fitted a Lorentzian (instead of the cut-off power law) to 
the same spectrum we could constrain neither the width 
nor centroid frequency. The fits were improved with the 
addition of the extra component, but were not sensitive 
to the exact form of that component. We fitted the spec- 
trum with a constant (instead of a cut-off power law or 
Lorentzian) and the fit also improved. None of the best-fit 
values of the peaked noise or VLFN were sensitive to the 
type of model that we used to describe the excess power; 
for simplicity, we therefore chose to use a constant. 

In the three highest intensity segments, the inclu- 
sion of a constant significantly improved the fits (deter- 
mined via an F-test for the inclusion of the extra compo- 



nent). Furthermore, in the segment with mean intensity 
3458 counts/s, fitting without a constant resulted only in 
determining an upper limit on the peaked noise of <3.1 % 
(Table ^). This upper limit is rather large, and a visual 
inspection of the power spectrum showed peaked noise 
to be present. Therefore, it seems justified to include an 
extra component in the fits, without needing (nor being 
able) to describe the actual form of the possible higher 
frequency feature. We note that: HFN was reported in 
Sco X-2, with cut-off frequencies of 66 Hz and 41 Hz in 
the NB and FB respectively (HK); also a QPO with a fre- 
quency of ~125 H z and FWHM of ~30 Hz was detected 
by Penninx et al. ( |1990D , in the FB of GX 17-1-2. 

The fact that we cannot determine the form of the 
higher frequency feature is not a major problem. The main 
result we obtained from the PC data is this: the peaked 
noise feature experienced a real weakening in the upper 
FB. This result is not sensitive to the presence or otherwise 
of a higher frequency feature. 



5. Discussion 

5.1. Observed properties of Sco X-2 

We were able to trace the power spectral features of Sco X- 
2 as a function of position in the Z track, in the energy 
range ~1-17 keV. At no time did we detect a narrow 
N/FBO as found in Sco X-1. Broad peaked noise was 
present at ~4-7 Hz, with FWHM ~6-12 Hz. We will 
subsequently refer to this feature as flaring branch noise 
(FBN). FBN was detected in the lower NB and reached 
a maximum at rank number 1.09 ('--^lO % of the way up 
the FB). At its maximum, the fractional rms was ^--^3 %. It 
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Table 4. Results from fits to the PC power spectra from the energy range 1.2-17.9 keV. 



Mean Intensity VLFN Peaked Noise Constant Reduced x /dof 



(cnts/s/4-dctcctors) Index rms (%) Frequency (Hz) FWHM (Hz) rms (%) 



Fits without an extra constant 
















3041±6 


-O.SitO.l 


l.liO.l 


5 4+0-9 


' ■'^-3.1 


2.5±0.4 




1.11 / 58 


3458±14 


-0.5±0.1 


1.2±0.1 


5.4" 


7.2" 


<3.1 




1.26 / 61 


4962±25 


-1.3±0.2 


2.0±0.3 


5.4" 


7.2" 


<1.8 




1.23 / 61 


5492±31 


-1.4±0.2 


9 7 + 0.5 
^' ' -0.4 


5.4" 


7.2" 


<1.3 




1.44 / 61 


Fits with an extra constant 














3041±6 


-1 0+"-^ 


l.OiO.l 


5-ot?:«o 


10.3+1° 


Q Q+0.4 

"^■"^-O.S 




1.04 / 57 


3458±14 


-1.7±0.3 


1 S + 

^■^'-o.s 




11.9tti 


Q 1+0.4 
•^■-^-0.3 


5.3e-6±1.0e-6 


0.58 / 57 


4962±25 


-1.5±0.2 


2.2±0.3 


5.9" 


11.9"" 


<1.6 


3.6o-6±l.le-6 


1.06 / 60 


5492±31 


-1.6±0.2 


o 1 +0.5 
■^■-^-0-4 


5.9" 


11.9" 


<1.2 


4.3c-6±0.9e-6 


1.10 / 60 



"Fixed parameter 



then became weaker and was undetectable by rank num- 
ber 1.5 (the middle of the FB; upper hmit 1.3 %). We 
found no significant change (more than a few Hz) in the 
centroid frequency or FWHM as the source moved up the 
FB. The FBN was stronger at higher energies. 

The behaviour we have observed is very similar to that 
reported by Ponman et at. ( |1988[ ). In addition, we have 
been able to unambiguously separate NB data from lower 
FB data. The increase in the fractional rms that accom- 
panied movement from the NB through to the lower FB, 
was also seen in the RXTE data of KK. We have not been 
able to confirm the decrease that they observed in the 
centroid frequency and FWHM, which accompanied that 
transition. 



5. 2. Comparison with other power spectral features 

The fast-time variability of Sco X-2 is different to that 
seen in the other Z sources. Most obvious is the low co- 
herence of FBN when compared with N/FBOs. The FBN 
has a relative width of ~1~2, compared with the N/FBO 
in Sco X-1 which has a relative width of less than '--^0.5 
(Dieters & van der Klis |200(j) . Also in Sco X-1, the fre- 
quency and FWHM of the NBO increase with increasing 
inferred mass accretion rate, and there is a sudden jump 
in frequency at the NB/FB vertex. In Sco X-2 on the other 
hand, there is a decrease in the FBN centroid frequency 
and FWHM with movement through the vertex (KK). In 
Sco X-1 the NBO strength increases with increasing rank 
number; the NBO becomes indistinguishable from the un- 
derlying noise (by ^10 % of the way up the FB) because 
the FWHM increases. In comparison, FBN is strongest 
at about 10 % of the way up the FB, and becomes unde- 
tectable because of a subsequent decrease in the fractional 
rms. 



Even given the marked differences, there are some sim- 
ilarities between FBN and the N/FBO seen in Sco X-1: 
the centroid frequencies are similar, both the N/FBO and 
FBN are hard, and there is no time lag between the soft 
and hard energies in either phenomenon. In Sco X-1, the 
NB/FB vertex is a critical position in the Z track. In 
Sco X-2 there are hints that this may also be true, al- 



though in a different sense (i.e. frequency jump for NBOs 
compared with a frequency fall for FBN). FBN and Z 
source N/FBOs may still be related, as already suggested 
by KK. 

In the scheme of the unified model, radiation pres- 
sure drives matter up out of the disc, forming an in- 



ner disc corona, resulting in a radial inflow (Lamb 1991 
Psahis et al. |1995D. When M rises to within -20 % of the 



Eddington limit, radial oscillations in the optical depth 
of the inflowing material produce a rocking of the X-ray 



spectrum; this gives rise to NBOs (Fortner et al. 1989 
Miller & Park |1995| ). This model may explain the evolu- 



tion of the N/FBO seen in Sco X-1, but it does not do well 
in explaining FBN. The FBN centroid frequency does not 
increase as predicted by the model. Furthermore, we have 
detected it —40 % of the way up the FB; oscillations are 
predicted to be suppressed at such high accretion rates. 

QPOs at super-Eddington accretion rates h ave been 
observed before: Kuulkers & van der Klis ( 1995b ) observed 
a -26 Hz QPO in the upper FB of Cyg X-2. The QPO 
was observed during an intensity dip in the FB, and it 
was stronger at lower energies. They suggested that it is 
produced by oscillations in an inner disc torus, and it is 
observable because Cyg X-2 is viewed at a high inclination 
(see Sect. 5.3). Sco X-2, on the other hand, is thought to 
be viewed at a low inclination (see Sect. |5.3| ); therefore it 
is unlikely that FBN is related to the 26 Hz QPO. 

KK pointed out that the FBN seen in Sco X-2 looks 
similar to peaked HEN sometimes seen in atoll sources. 
We observed a decrease in the FBN fractional rms as the 
source moved up the NB, which is opposite to what would 
be expected if the FBN was infact atoll source HFN. The 
low Nyquist frequency is no problem here as the break fre- 
quency of HFN decreases with decreasing inferred mass- 



accretion rate (Wijnands & van der Klis 1999). The phe- 
nomenology of FBN is thus significantly different to atoll 
source HFN; therefore we propose that FBN is unrelated 
to atoll source HFN. 
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5.3. Two types of Z sources: a clue to the origin of FBN? 

In order to understand why Sco X-2 behaves so differently, 
it is worthwhile to investigate the general relationships 
that exist between the Z sources. On the basis of their Z 
track morphology and fast-time variability, they can be di- 
vided into two groups: Cyg X-2, GX 5-1, GX MO+0 (Cyg- 
hke objects); and Sco X-1, Sco X -2, GX 17+2 (Sco-like 
objects) (e.g. Kuulkers et al. 1997 ). The Cyg-like sources 
have much larger HBs and smaller FBs than the Sco-like 
sources. The HBs in the Sco-like objects are slanted, while 
the HBs in the Cyg-like sources are horizontal. Most im- 
portantly, in the Cyg-like objects the intensity is seen to 
fall in the middle and upper FB, while the Sco-like objects 
really do flare in the FB. 

At high accretion rates an inner disc torus may be 
produced, which in high inclination sources can partly ob- 
scure the X-ray emitting region. The Cyg-like objects have 
thus been interpreted as being viewed at higher inclina- 
tions, and t he Sc o-like objects at lower inclinations (Ku- 
ulkers et al. 1997 ). In support of this, optical observations 
place Cyg X -2 at an inclination of 62.5±4° (Orosz fc Ku - 
ulkers |l999|) , Sco X-1 at 15°-40° (Crampton et al. |l976|) , 
and also suggest that Sc o X-2 is viewed at a l ow in clina- 
tion (Wachter & Margon |l996l ) . Psaltis et al. ( |l995|) have 
modelled the X-ray spectra of Z sources. They suggested 
that the Sco-like objects have weaker magnetic fields than 
the Cyg-like objects. 

The most obvious difference between Sco X-2 and the 
other Z sources is the absence of an HB, so it is tempting to 
think that either: Sco X-2 has a weaker magnetic field than 
Sco X-1, while also (perhaps incidentally) being viewed at 
a low inclination; or simply that Sco X-2 never reaches a 
low enough mass- accretion rate to trace out an HB. 

A close investigation of the timing properties of Cyg X- 
2 may provide hints as to the origin of FBN. In Cyg X- 
2 long-term variations in the mean X-ray intensity have 
been discovered (e.g. Paul et al. 2000| ). Ginga observations 
of Cyg X-2 showed that the Z track was more Sco-like 
at medium level intensities, and more Cyg-like at high 



level intensities (Wijnands et al. 1997). RXTE observa- 



tions of Cyg X-2 at low intensities revealed no NBO (up- 
per limit ~1 %) in the NB; however a peaked noise fea- 
ture was detected in the lower FB, with a frequency and 
fractional rms of 6-7 Hz and ^^3 % respectively (Kuulk- 



ers et al. 1999 ). This is (perhaps only superficially) similar 
to the FBN seen in Sco X-2. 

It is vital to have more observations of Cyg X-2 at low 
intensities so that we can better understand the proper- 
ties of the peaked noise feature. Only then can we make 
any conclusions about the relationship between FBN and 
the Cyg X-2 peaked noise. But if the two features are the 
same phenomenon, we must ask: what are the physical 
characteristics of Cyg X-2 when it is at low level inten- 
sities? The low level may be due to obscuration by the 
outer accretion disc, possibly caused by a precessing disc 



(e.g. Kuulkers et al. |996|). Wijnands & van der Khs ( |200lD 
have argued against the precessing disc model on the basis 
of differing kHz QPO properties at different intensity lev- 



els. Kuulkers et al. (1999) showed that the relative VLFN 



properties at different intensity levels were also inconsis- 
tent with obscuration. Instead, they preferred a model in 
which the overall intensity level is sensitive to the proper- 
ties of the inner mass flow. Whatever the case, if FBN and 
Cyg X-2 peaked noise are produced by the same process, 
then that process will be largely insensitive to orbital in- 
clination. In this way, Sco X-2 may prove to be a useful 
tool in discerning between models trying to explain the 
behaviour of Cyg X-2, and ultimately the other Z sources. 

6. Conclusions 

We have carried out the first detailed study of the X-ray 
fast-time variability of the Z source Sco X-2, as a func- 
tion of position in the Z track. We found broad peaked 
noise with centroid frequency and FWHM in the range 
~4-7 Hz and ~6-12 Hz respectively. We find it difficult 
to attribute FBN to the same processes which are able to 
explain typical Z source N/FBOs. We also conclude that 
it is not a manifestation of atoll source HFN. FBN resem- 
bles the FB peaked noise seen in the FB of Cyg X-2 at low 
overall intensities. If these are later confirmed to be the 
same phenomena, then the process responsible for FBN 
must be largely insensitive to orbital inclination. 
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Appendix: Calculation of deadtime-affected virhite- 
noise levels 

In the absence of deadtime, the white-noise level in a 
Leahy normalised power spectrum, due to Poisson noise 
in the data, is equal to 2 (Leahy et al. 1983|) . The effect of 
deadtime is to decrease the variance of the data and thus 
lower the white-noise level. The deadtime-affected white- 
noise level for the sum channel P^s, can be calculated for 



each frequency v in the power spectrum 



< P,.s >= 2(1 -^sTdcad) X 

Ms'^dcad \ / ^dcad 



1 



1 /^s'^dcad 



sm 



\ 2z^Nyquist 



(A.1) 
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where /is is the observed counting rate per detector, Tres is References 



the time resolution of the lightcurve, Tdoad is the deadtime 
per event, and ^isryquist is the Nyquist frequency of each 
transform (Weisskopf 1985| ). 

The deadtime-affected white-noise level for any partic- 
ular energy band P^h, can be calculated as a function of 
Pys- For an arbitrary time interval, the ratio of the total 
number of counts k in the particular energy band, to the 
total number of counts n in the sum channel is denoted 
by r. The mean variance of the data is denoted as and 
(t2 for the sum channel and the particular energy band re- 



spectively. These are related by (Mitsuda & Dotani [1989D 



2 2 



+ k{l-r) 



(A.2) 



Each transform is calculated from a time series Xk 
(k = 0, ...,Af — 1) consisting of N arbitrary time inter- 
vals. By employing Parseval's theorem, we can relate the 
total variance Var{xk) (or the mean variance cr^) of the 
time series, and the Fourier amplitude a 



Var{xk) = Na"^ 



N_l 

N ' 



(A.3) 



Here, we have considered only situations where the depen- 
dence of Fourier amplitude on frequency is negligible (i.e. 
white-noise). The factor of N-1 originates from the fact 
that the sum used in the expression for the total variance 
does not include a term for a frequency of zero. 

The mean variance and Leahy power P can then be 
related via 



a'^P- 



iVph fN-l 



iV2 



(A.4) 



A^ph is equal to Nn for the sum channel or Nk for the 
particular energy band. We can express Pb as a function 
of Ps by substituting Eq. into Eq. ^3. 



Ph- 



k f N -I 



N 



^'n2 2 



k'^ n f N -I 



N 



+ k{l-r) (A.5) 



or 



Ph = 2 



N 



^A- 1 

which in the case of large A^ reduces to 



Pb = 2 - - (2 - Ps) 
n 



Ps 



(A.6) 



(A.7) 
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